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0\ ■ ABSTRACT 



I review the results from lattice gauge theory for the masses of exotic hybrid mesons. 



1 Introduction 



The quark model predicts that the charge conjugation (C) and parity (P) of a meson 
with spin S and orbital angular momentum L is 

P=(-1) L+1 C=(-l) L+s (1) 

States with quantum numbers not produced by eq. [I], such as 

^ c = 1^0 + -,2+-,0- (2) 

are known as exotics 0). Exotic states are allowed by QCD. 

There are a number of different possibilities for the structure of an exotic 
state. An exotic signal could be: a hybrid meson, which is quark and anti-quark and 
excited glue, or bound state of two quarks and two anti-quarks (QQQQ). The two 
most popular guesses for the structure of the (QQQQ) state are either a molecule of 
two mesons or diquark anti-diquark bound state. In this paper I review the latest 
lattice results for the masses of exotic hybrid mesons, concentrating on the 1~ + 
state, obtained from lattice QCD. 



2 Lattice simulations of exotic mesons 



Many numerical predictions of QCD can be determined from the path integral 

c{t) ~ J dU J dip J dip J2O(0,0)O(x,ty e - SF - SG (3) 

where Sf is the fermion action (some lattice discretization of the Dirac action) and 
Sq is the pure gauge action. The path integral in eq. ^] is put on the computer 
using a clever finite difference formalism @), due to Wilson, that maintains gauge 
invariance. The physical picture for eq. [I] is that a hadron is created at time 0, from 
where it propagates to the time t, where it is destroyed. The fermion integration 
can be done exactly in eq. ||] to produce the fermion determinant. Simulations that 
include the effect of the determinant are very expensive computationally, so typically 
it is not included in the simulation (the quenched approximation). However there 
has been some recent work that includes the effect of the determinant 0) on the 
light exotic spectrum. 

The standard interpolating operator for the pion, which can be used in 

eq. |, is 

O n (x,t) =fet) 75 ^(x,t) (4) 

which has the correct J PC = h quantum numbers. One possible interpolating 
operator ED for an exotic 1~ + particle is 

0-L-+(x,t) = ip(x,t)"/jF i:j (x,t)ip(x,t) (5) 

where F is the QCD field strength tensor. It is essential to use operators that 
have some spatial separation between the quarks in the meson to get a good signal. 
Recently the MILC collaboration has attempted to measure the "wave function" of 
the 1 h hybrid meson in coulomb gauge HP. Unfortunately the operator used did 
not have the correct charge conjugation quantum number, so the published wave 
function §P is incorrect. 

In this formalism a gauge invariant interpolating operator, for any possible 
exotic hybrid particle or four particle state can be constructed. The dynamics then 
determines whether the resulting state has a narrow decay width, which can be 
detected experimentally. In the large N c (number of colours) limit E> H) both exotic 
hybrid mesons and non-exotic mesons have widths vanishingly small compared to 
their masses. 

The data from the simulation is extracted using a fit model @) : 

c(t) = a exp(—m t) + aiexp(-mit) + ■ ■ ■ (6) 



where itlq (mi) is the ground (first excited) state mass and the dots represent higher 
excitations. Although in principle excited state masses can be extracted from a 
multiple exponential fit, in practice this is numerically non-trivial because of the 
noise in the data from the simulation. Simulations that involve the calculation of 
the properties of exotic hybrid mesons are harder than those that concentrate on 
the non-exotic hadrons, because the signal to noise ratio is worse for exotic mesons 
than for QQ mesons. 

In an individual lattice simulation there are errors from the finite size of 
the lattice spacing and the finite lattice volume. State of the art lattice simulations 
in the quenched theory, run at a number of different lattice spacings and physical 
volumes and extrapolate the results to the continuum and infinite volume. For 
the exotic mesons, this is done for heavy quarks (see section |J), but as yet, the 
continuum extrapolation has not been done for light exotic mesons. 

How successful is lattice QCD in practice? One way to tell is to compare 
the predictions for the masses from lattice QCD for well known particles (proton, p, 
etc.) with experiment. The most accurate quenched calculation to date has recently 
been completed by the CP-PACS collaboration 0). From the masses of 11 light 
hadrons, they conclude that the quenched approximation disagrees with experiment 
by at most 11%. The comparison of results from simulations that include dynamical 
fermions with experiment is less clear, because of their high computational cost (see 
Kenway ^ for a review of the latest results). 

The results from lattice QCD also provide insight into the underlying dy- 
namics of light hadrons. Lattice QCD simulations can test the various assumptions 
made in models of the QCD dynamics. For example there are a number of models 
of exotic states based on the idea of a bound diquark anti-diquark pair 0) . A crit- 
ical assumption in diquark models is that two quarks actually do cluster to form a 
diquark. This assumption has recently been tested in a lattice gauge theory simu- 
lation by the Bielefeld group EHD , where they found no deeply bound diquark state 
in Landau gauge. 

3 Results for light exotic mesons 

In the last year the MILC collaboration have repeated ID their initial simulations H) 
using the clover fermion action. The clover action is "closer" to the continuum 
than the Wilson fermion action, because it has the leading order lattice spacing 
terms removed. There are also new results for the hybrid masses from the SESAM 
collaboration (reported by Lacock and Schilling) 0) , that include some effects from 



Table 1: Masses of the light 1 + hybrid from lattice gauge theory. 



Simulation 


Group 


mass GeV 


A 
B 
C 
D 


UKQCD ED 
MILC i D 
MILC i § 
Lacock and Schilling 0) 


2.0(2) 
2.0(1) ±sys 
2.1(1) ±sys 

1.9(2) 



Table 2: Parameters of light exotic meson simulations, 



Simulation 


Action 


fermions 


length (fm) 


a" 1 GeV 


MpsjMy 


A 


clover 


quenched 


1.6 


2.0 


0.76 


B 


Wilson 


quenched 


2.3 


2.8 


0.96,0.93,0.88,0.77 


C 


clover 


quenched 


2.3 


2.8 


0.94,0.90,0.72 


D 


Wilson 


dynamical 


1.4 


2.3 


0.83,0.81,0.76,0.69 



dynamical sea fermions. 

The results for the mass of the 1~ + exotic state are summarised in table [I]. 
All the results from the various simulations are essentially consistent, with the mass 
of the 1 h state around 2 GeV. In table ^| we show the physical parameters for each 
of the simulations. The interpolating operators used to create the exotic meson 
states in the MILC calculations © are different to those used in the UKQCD 0) 
and SESAM simulations 

The observation that the results for the mass of the 1 h hybrid meson 
are consistent for very different simulations gives us confidence in the final re- 
sult. Although I would prefer to see simulation results at lighter quark masses. 
Simulations at the point where the ratio of the pseudoscalar mass to vector mass 
(Mps/My ~ 0.5) are possible with current algorithms and computers 0). 



4 Results for heavy exotic mesons 

There has been a lot of work on calculating the spectroscopy of cc and bb mesons from 



lattice gauge theory (see Davies for a review) . The main technical complication 
in heavy quark simulations is that the lattice spacing of current simulations is not 
smaller than the heavy quark mass. So various effective field theory Lagrangian 
approximations to QCD are simulated. 

The NRQCD (nonrelativistic QCD) Lagrangian is one such effective field 
theory approximation to QCD, with the expansion parameter equal to the veloc- 



Table 3: Mass splitting between the bb 1 + hybrid and the bb IS state. 



Group 


comments 


mass GeV 


UKQCD 16 P 
CP-PACS Pi) 
Juge et al. ^) 


0(Mv 4 ) errors 
Asymmetric, 0(Mv A ) errors 
Asymmetric, 0(Mv A ) errors 


1.68(10) 
1.542(8) 
1.49(2)(5) 



ity squared. NRQCD has been particularly successful in simulating the Upsilon 
spectrum 0), but is less well converged for charmomium, (particularly for spin 
splittings), because the charm quarks move with higher velocity 0). The NRQCD 
Langrangian correct up to 0(Mv 2 ) is 



C NRQCD = ^_ 



A 2 c A 4 Cl aB 



(7) 



2M 8M 2 2M 

where c and C\ are coefficients obtained by a perturbative matching procedure to 
QCD. In table [3] the results of all the recent NRQCD simulations of the bbg hybrids 
in the quenched approximation are compiled. In the hybrid meson simulations no 
spin terms are included in the Lagrangian (ci = in eq. |7|), so the 1~ + , H , and 



TJ) 



2 H states are degenerate. Both the results from the CP-PACS collaboration 
and from Juge, Kuti and Morningstar HID were shown to be independent of the 
lattice spacing and lattice volume. For example, the CP-PACS collaboration 0) 
found that the masses of the hybrids were independent of the box size above 1.2 fm. 

The "asymmetric" comment in table |3| refers to the technique of treating 
space and time asymmetrically. A smaller lattice spacing was used in the time 
direction than in the space direction, which allowed the signal to be seen for further, 
for a given spatial volume. This technique has helped to reduce the errors. The 
practicalities of this idea were demonstrated by Morningstar and Peardon for 
the glueball spectrum. 

In table |] the results of the mass splitting of the 1 _+ states and the IS state 
are shown in the charmonium system. The MILC collaboration used the standard 
Wilson and clover actions to simulate the charm quark in their simulations of heavy 
exotic mesons, as previous work has shown this to be reliable PD 

The first results for heavy exotic EZP hybrids were done in the adiabatic 
surfaces approach, where the effect of the excited glue is subsumed in a potential 
measured on the lattice (see for a review). Juge, Kuti and Morningstar li^J 

have completed a systematic study of these potentials. The NRQCD approach is a 
more accurate approximation to QCD than the adiabatic potential technique; how- 



ever Juge, Kuti and Morningstar ED found that the adiabatic potential approach 



Table 4: Mass splitting between the cc 1 + hybrid and the cc IS state. 



Group 


comments 


mass MeV 


MILC S § 

MILC I) 
CP-PACS El) 


Wilson action 
clover action 
NRQCD 0(Mv 4 ) errors, Asymmetric 


1340±?° + sys 
1220t$j + sys 
1323(13) 



reproduced the level splittings from NRQCD up to 10 %. A preliminary result for 
the calculation of the adiabatic surfaces with the effects of dynamical fermions in- 
cluded has been reported by Bali 0) and collaborators. No dramatic differences 
between the quenched theory result were observed. 



5 Conclusions 

All the lattice simulations agree that the light 1~ + state has a mass of 2.0(2)GeV. 
The first simulation that included the effects of dynamical fermions has not changed 
the result. 

The experimental results for light exotics are reviewed by S.U. Chung 0), 
so I just briefly compare the lattice results to experiment. There is an experimental 
signal for a T + state at 1.4 GeV with a decay into r/ir EHI) from E852, Crystal 
barrel, VES, and KEK. It is surprising that this state is only seen in the f]ir channel, 
as this decay is theoretically suppressed relative to other decays The E852 

collaboration have also reported H3) a signal for 1 _+ state decaying into p-n with a 
mass of around 1.6 GeV. The decay width is in reasonable agreement with theoretical 
calculations EH) . 

Clearly the agreement between the possible experimental signals for the 
1~ + states and the lattice results is very poor. The errors on the lattice results for 
1~ + states are large relative to the errors on QQ states. To quantify the disagreement 
between experiment and lattice results the systematic errors on the lattice simulation 
results should be reduced. In particular the masses of the quarks used in the lattice 
simulation should be reduced. It is possible that the states seen experimentally are 
really QQQQ states, in which case the operators used in the lattice simulations 
would not couple strongly to them. 

Although the adiabatic lattice potential approach is not expected to be a 
good description of the physics of light hybrids, we note that the results of Juge, 
Kuti and Morningstar , show that the splitting between the ground and first 
excited state is about 200 MeV, in broad agreement with the experimental results 



of the E852 collaboration. Although no insight is gained about the different decay 
widths. 

To definitely identify an exotic hybrid meson requires both the calculation 
of the mass as well as the decay widths. There has been very little work on hadronic 
decays on the lattice. The most obvious hadronic process to study using lattice 
gauge theory is the p — > tttt decay, however there have only been a few attempts to 



calculate the g pn7r coupling ED The GF11 lattice group has recently computed the 
decay widths for the decay of the 0++ glueball to two pseudoscalars H. 

The MILC collaboration @) have started to investigate the mixing between 
the operator in eq. § and the operator (tt <8> ai) eq. [|. 

?W°?757^ 6 (8) 

which has the quantum numbers This type of correlator would naively be 

expected to yield the decay width of the 1~ + state to p, and a\. Unfortunately 
the analysis of Maiani and Testa EH) shows that the matrix element required in 
the computation of the decay width is hidden beneath an unphysical term that 
increases exponentially with time. The origin of the unphysical term comes from the 
requirement that both final mesons should be onshell and is deeply related to theory 
being defined in Euclidean space (required for us to have a well defined theory to 
simulate). Some information may be extracted for onshell processes EH El), using 
the methods proposed by Michael E3). 
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